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ABSTRACT 

Compilations of abundances of Strontium and Barium in extremely metal- 
poor stars show that an apparent cutoff is observed for [Sr/Ba] at [Fe/H]<-3.6 
and large fluctuations for [Fe/H]>-3.6 with a clear upper bound depending on 
metallicity. We study the factors that place upper limits on the logarithmic ratio 
[Sr/Ba]. A model is developed in which the collapses of type II supernovae are 
found to reproduce many of the features seen in the data. This model is consistent 
with galactic chemical evolution constraints of light-element enrichment in metal- 
poor stars. Effects of turbulence in an explosive site have also been simulated, 
and are found to be important in explaining the large scatter observed in the 
[Sr/Ba] data. 

Subject headings: stars: Population II — nuclear reactions, nucleosynthesis, 
abundances — black hole physics — Galaxy: evolution — equation of state 
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Introduction 



Nuclei heavier than iron and nickel are mostly made either in the s-process (s. 



neutron-capture proce ss) or r-process (rapid neutron-capture process) f lB2FH 



1957 



ow 



Wallerstein et al 



19971 ). which is thought to produce half the nuclides heavier than Fe and 



nearly all actinides observed in m ost metal-poor stars- 



can o ccur in first generation stars ( iGoriely fc Siess 



2001 



Although effects o 



Suda et al. 



2004; 



the s-process 



Nishimura et al. 



20091 ). the s-process appears not to be responsible for the enrichment of the interstellar 
medium with the elements produced in the early life of the galaxy because it requires longer 
timescales to populate the ISM with its products. The r-process, however, produces its 
nuclei earlier in the galactic history. Because massive stars produce most of the required 
conditions for the r-process, th ey have long been s upposed to synthesize its nuchdes. 



Howe ver, detailed calculations ( iFischer et al. 



2010 



Hudepohl et al.l 



2010 



Goriely et al. 



20111) have suggested that massive stars may not fulfill all of the r-process conditions. We 
present data that sheds important new light on this question. 

Stars with masses greater than 8 Mq are thought to become core-collapse supernovae. 
Those less than 25 Mq will form neutron stars, emitting neutrinos and driving some of 
their newly synthesized nuclei into the interstellar medium. Those more massive than 
about 40 Mq are thought to collapse directly to black holes, contributing nothing to the 
interstellar medium. Stars with 25< M < 40 Mq are thought to collapse to neutron stars 



then, as some of t 



2003 



le expelled material falls back onto the neutron star, to black holes ( Fryer 



Heger et al 



20031 ). producing "accretion-induced black holes" and, assuming some 



material can escape, would be expected to contribute some newly synthesized nuclides to 
the interstellar medium. 

We report on some striking astronomical data from extremely metal-poor (EMP) stars, 
those having [Fe/H]<-2.5. Their metals are thought to have been produced in the first 
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generation massive stars that were the progenitors of the ones currently observed. Because 
these stars have very low abundances of anything except primordial hydrogen and helium, 
their abundances of heavy nuclei are presumed to arise primarily from the r-process. We 
concentrate on elements made in the r-process and that are easy to detect in these stars: 
strontium and barium. 



Previous models ( IQian fc Wasserburg 



2008 



Wanajo et al. 



2011 



Travaglio et al. 



20041 ) 



have described separate r-process sites as the cause of an enrichment of light neutron-capture 
r-process elements. We contrast our model with the "weak" r-process and the light element 



primary process (LEP 



elements (IHonda et al. 



) by describing an enrichment in the lighter neutron-capture 



20061 ) that is a natural consequence of a single site but with differing 



levels of enrichment of light and heavier neutron-capture elements. Our model produces 
the r-process from an explosive event (e.g., a type II supernova), while the r-process ejecta 
is prohibited in some fraction of those events from escaping due to a subsequent collapse 
of the proto-neutron star to a black hole. The collapse time in this event relative to the 
r-process ejection time will determine how much r-process material is ejected, as well as its 
composition. 



Observational Data 



20081 ). The sample includes 



Data were taken from the SAGA Database (iSuda et al. 
data from 186 papers and 4270 stars (1488 unique stars) dealing with observations of 
metal-poor halo stars published prior to 2012. The number of stars for which [Fe/H]<-2.5 
is 428, reducing to 361 stars if carbon-enhanced EMP stars are excluded. This sample 
could suffer errors due to systematic offsets between different studies; a comparison of 
results for the same stars but from different studies shows that these can sometimes be 
as large as a factor of 3. We thus selected the results of Sr and Ba from observations in 
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which both elements are reported in the same paper - 260 total stars. The abundances in 
this database were all obtained assuming local thermodynamic equilibrium (LTE). The 
scatter found in Sr/Ba as well as in Sr/Fe and Ba/Fe in that sample in the metallicity 
range of -3.6<[Fe/H]<-3.0 is observed to be a factor of 100 or more, much larger than the 
observational errors. 

Figure [1] shows the astronomical data for [Sr/Fe], [Ba/Fe], and [Sr/Ba] from the data 
set described above, showing a wide range of [Sr/Ba], [Sr/Fe], and [Ba/Fe] values at each 
value of [Fe/H]. Our interpretation of these features will be given in subsequent sections; 
the curves in Figure [T] that overlay the data, reflect these efforts. 

We find three remarkable features in the [Sr/Ba] distribution. 



1. Almost all stars appear in the range of the functional [Sr/Ba]< —1.5— [Fe/H]. There 
are two exceptions in [Fe/H]< -2.0. One of them is BD+80"245 ([Fe/H]= -2.09, 
[Sr/Ba]= +1.04), which is a well known " a-deficient" star (e.g., [Mg/Fe]= -0.22) 



having exceptionally low abundances o: 
[Ba/Fe] = —1.89 which wa s stud ied by 



neutron-capture elements ([Sr/Fe]= —0.85 and 



Ivans et al. 



studied by 



Barklem et al 



( 120031 )). The other is HE0029-1839 



( 120051 ). for which little information exists. 



2. There are only a few stars with [Sr/Ba]<-0.5. We present calculations that show that 
this range is not expected from the r-process or from tr-process, which we discuss in 
the next section. The exceptions follows: 

• CS 30322-023 ([Fe/H]=-3.25, [Sr/Ba]=-1.10): a moderatel y carbon-enhan ced 



star ([C/Fe 



Aoki et al. 



0.6) with enhancement of nitrogen and Ba ( iMasseron et al. 



2006 



20071 ). 



• CS 29493-090 ([Fe/H]=-2.82, [Sr/Ba]= -1.41): a moderately carbon-enhanced 
star ([C/Fe]=+0.74) with excess of Ba ^Barklem et al.l 120051 ). 
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• CS 22946-011 ([Fe/H]=-2.59, [Sr/Ba]=-0.99) and CS 22941-005 ([Fe/H]=-2.43, 
[Sr/Ba]=-0.7 7): binaries with excesse s of Ba ([Ba/Fe]=+1.26 and +0.34, 



respectively: 



Preston fc Snedenl toom . CS 22950-173 ([Fe/H]=-2.5, [Sr/Ba] 



0.72) which also 
([Ba/Fe]=-0.04: 



belongs to be a b i nary s ystem, though Ba is not enhanced 



Preston fc Snedenl (120001 )). The CH feature would not be 



detected due to high temperature (~6800K). 

HE0305-4520 ([Fe/H]=-2.91, [Sr/Ba]=-1.25): Ba is slightly enhanced 



([Ba/Fe]=+0.59), w: 



Barklem et al. 



lile an excess of carbon is not clear ([C/Fe]=+0.33: 



mom ). 



Hence, except for HE0305-4520, all the above stars are (candidate) carbon-enhanced 
stars, which can be affected by an s-process contribution through mass transfer in 
binary systems. HE0305-4520 could also be slightly affected by s-process. In the 
objects with low Sr abundances, Sr/Ba could be driven to very low values even by a 
small contamination of Ba from the s-process. 

There is no object for which s-process contamination is not suspected in the region 
[Sr/Ba] < —0.5 given the typical errors (~0.2 dex or larger) in [Sr/Ba]. 

3. A cutoff in [Sr/Ba] data at [Fe/H]~ —3.6, which will be addressed below, is 
found. There are only four stars (CS 22172-002, CD-38°245, CS 22885-296, and 
CS 30336-049) in [Fe/H]< —3.6 for which both of the Sr and Ba abundances are 
determined, and all of them have low [Sr/Ba] (< 0.05). There are roughly three 
times as many stars with [Sr/Ba]> 0.05 than with [Sr/Ba]< 0.05 in the region 
—3.6 <[Fe/H]< —3.0. If the same distribution is assumed for —4.0 <[Fe/H]< —3.6, 
one would expect 12 - 15 stars with [Sr/Ba]>0.05, instead of the zero stars that are 
observed, from the four stars with [Sr/Ba] >0. 05. 

To confirm this cutoff, we inspected the data with [Fe/H]< —3.6 in more detail. We 
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found 18 objects in this range in our sample, but both Sr and Ba abundances were 
determined only for seven objects, among which three are carbon-enhanced objects. 

Among the remaining 11 objects for which [Fe/H]< —3.6, five stars are carbon- 
enhanced, and four have very low [Sr/Fe] or low upper limits, suggesting that Sr/Ba 
ratios are also quite low. There is no signature to have high Sr/Ba for other two 
objects although the upper limits for Ba abundances are still weak. 

We note that no new stars with both Sr and Ba measurements for [Fe/H]< —3.6 were 
re ported by th e recen t paper to study chemical composition of most metal-poor stars 



by 



Yong et al 



mm . 



3. Galactic Chemical Evolution Interpretation and Turbulence 



The most daunting challenges for any attempt to explain the data distribution for the 
EMP stars are reproduction of the sharp cutoff in [Sr/Ba] at [Fe/H]=-3.6, the explanation 
of the upper and lower limits as a function of metallicity in that distribution, and the wide 
dispersion in [Sr/Ba]. 

A red uction in the Sr a nd Ba production in a SN collapse sc enario was studied 



previously (IBoyd et al. 



2012), using the formalism and results from 



Woosley et al. 



(ll99J) 



There, the r-process occurs in upper-most 15 shells that appear just at the surface of the 
nascent neutron star. The mass, the r mody namic trajectories, and initial composition of each 
shell are described in 



Woosley et al 



( 119941 ) . Those shells were examined in the framework of 
the ejection time, the black-hole collapse time, and the resultant nucleosynthesis, assuming 
that the ejection is halted by the collapse to a black hole. We assume that any collapse time 
occurring prior to the shell ejection time in this model will prevent that shell from being 
ejected. At present, we extend that model to study the relationship between the stellar 
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metallicity, the progenitor mass, the collapse time, a nd elemental yields. This model is 



19951). 



coupled to a galactic chemical evolution (GCE) code ( iTimmes. Woosley. fc Weaver! 
The net result is that a more massive progenitor will generally collapse earlier, resulting in 
ejecta enriched in the lighter neutron-capture elements (A<130) but a lower yield of the 
heavier elements. Th is GCE formulation will be described in detail in a subsequent paper 



( iFamiano et al.l 



2OI2I ). Here, we present the general principles. 



Previous results have indicated that a truncation in the r-process in type II supernovae 
due to a collapse of the proto-neutron star to a black hole may be responsible for a relative 
enrichment of light neutron-capture elements in som e metal-poor halo stars produced by 



a single-site "truncated r-process" ( iBoyd et al. 



2OI2I ) (where single-site is defined to be 



a single episode of star formation after one r-process event). This "tr-process" could be 
caused by either a dynamic collapse of material below the black hole mass cut or by a 
change in the neutrino luminosity resulting in an increased electron fraction Y^. in the later 
stages of the r-process, which would result in a reduced synthesis of the r-process nuclei 
with masses in excess of A~130. 



20121 ) in which the isotopic 



We adopt a model based on prior calculations ( iBoyd et al.l 
yields in a tr-process are directly related to the collapse time in an accretion-induced black 
hole. In order to produce the chemical evolution results shown, we began with the Ba 



yields of 



Cescutti et al 



( I2OO6I ) and performed several calculations iteratively to arrive at 
the yields shown in Figure |2l These assumed yields reproduce the GCE results shown by 
the red lines in Figure [TJ For stars with 10 Mm <M< 30 Mq, the yields are roughtly equal 



to those of Model IMax in 



Cescutti et al. 



in the same reference and less than those of 



(120061) but slight 



Raiteri et al. 



y less than those of model 2Max 



(119991 ) by about a factor of four 



(for 1 0-11 Mq SNIIs). For M<10 Mq yields are estimated from Figure 3 of 



Cescutti et al. 



( I2OO6I ). The Sr yields (Figure [2]) were chosen to reproduce the [Sr/Fe] and [Sr/Ba] GCE 
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results i n Figure [H The 



those of 



Cescutti et al. 



only significant difference between the Ba yields in this paper and 



(120061 ) corresponds to 30 Mq stars, which are roughly 30 times 
higher than those of model IMax. This produces Ba early in the GCE model, though the 
contribution is insignificant at late times. 

A power law IMF was used with a Salpeter exponent of -1.31. The black hole collapse 
time was cor nputed based on progen itor mass and metallicity and the nuclear equation of 



state (EOS) flO'Connor fc Ott 



20101 ). This time was then assumed to be the truncation 



time in a dynamic tr-process. For later truncation times, mass shells closer to the 
proto-neutron star core, blown off later in the explosion model, are assumed to be ejected 
either by achieving a velocity in excess of the escape velocity or perhaps by being ejected 
in a jet. Later collapse times correspond to larger production yields of heavier r-process 
elements. 

In the r-process, most of the Sr is produced in the early-stage ejectt in mass shells 
ejected within the first four seconds post-bounce, while most of the Ba is produced in 
the later-stage ejecta of a SNII in mass shells ejected from between four and 18 seconds 
post-bounce. Thus, the Sr yields in a tr-process are fairly constant with collapse time and 
fairly close to their complete r-process values except for very short collapse times for which 
no r-process material is ejected. As only fairly late collapse times will allow for a significant 
amount of Ba ejection in the tr-process, [Sr/Ba] approaches zero in the later-stage ejecta. 



2010l ) were used to 



Results from the spherical collapse code GRID ( 10 'Connor fc Ott 
study the relationship between the progenitor star's mass and metallicity and the Sr and 
Ba yields in massive stars. It is noted that these yields are minimum yields, since in the 
non-rotating spherical models employed here, the collapse times are minimum collapse 
times. Longer collapse times may result from effects such as asymmetric explosions, 
rotations, and neutrino heating induced by these effects. However, since Sr is produced in 
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shells that are ejected very early on, increasing the collapse time will have little effect on 
the Sr production, while significantly increasing the Ba production. 

Ejection of supernova materials into the interstellar medium resulting from turbulence 
in the stellar interior may also affect the observed abundances. Although it is difficult to 
simulate such effects in a one-dimensional model, we have performed a simplified analysis 
that at least suggests the sort of effects that might be observed, and probably to constrain 
the magnitude of the predictions. This analysis took each of the 15 individual mass shells 
that were used in our analysis and that could be ejected from the surface of the neutron 
star. We then assumed that turbulent ejection might allow a single individual shell to be 
ejected from the star, while the others collapsed onto the proto-neutron star surface. As 
noted above, the outer shell would be rich in Sr and would have very little Ba, while the 
inner shells would have much less Sr and would be much richer in Ba. Ejection of only the 
top or bottom shells should provide limits on the ability of turbulence to explain the range 
of [Sr/Fe], [Ba/Fe], and [Sr/Ba] observed and might also explain the dispersion in the data 
that mixing might produce. The [Sr/Ba] results from sites ejecting only these shells are 
shown in the GCE model in Figure [T]i. 



Results 



Assuming a tr-process model resulting in a reduced Sr and Ba production, the 
single-site [Sr/Ba] values as a function of metallicity [Fe/H] are sho wn in Figure [1] with a 



collapse model utilizing the LS220 EOS OLattimer fc Swesty 



2004j ). which is one of the 



most widely used EOSs in numerical simulations of supernova explosions. Several other 
EOSs have been examined, and it has been found that a softer EOS results in progressively 
larger [Sr/Ba] values at lower metallicity (corr esponding to earlier times in the galactic 



evolution produced by stars with M> 20 Mq) (IFamiano et al. 



2OI2I) 
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In Figure [H two GCE results are presented. One is the metallicity relationship from a 
model based on the input yields from Figure [2] showing [Sr/Fe], [Ba/Fe], and [Sr/Ba] ratios. 
It can be seen that those fall within the extremes of the observational data points. The 
other result is a metallicity relationship corresponding to a GCE model in which all stars 
with M>2OM0 ultimately collapse to black holes. 

In this model, an earlier collapse means that shells closer to the neutron star surface 
(which produce more Ba) will not be ejected, and the overall Ba ejection (summed over all 
shells) is reduced. Since Sr is mostly produced in outer shells, it is not reduced as much. 
The net result is that the calculated lines in the figure will correspond to the minimum in 
the [Sr/Fe] and [Ba/Fe] values and a maximum in the [Sr/Ba] values. Collapse times range 
from about 0.25s for 40 Mq stars with a soft EOS to over 3.5s for 20 Mq stars with a stiff 
EOS. 

The results of these calculations can be compared to the astronomical data of [Sr/Ba] 
discussed here. These observations show a sharp maximum in the [Sr/Ba] value for a given 
metallicity. This could be explained by the tr-process limiting the Sr or Ba production as 
a function of stellar age or mass. We also note a similarly sharp lower limit in the [Ba/Fe] 
values at [Ba/Fe]<0 corresponding to the sharp [Sr/Ba] upper limit. However, the lower 
limit in [Sr/Fe] is less distinct. The curves shown in Figure [Hi represent attempts to study 
the effects of turbulence by assuming that individual mass shells were ejected into the ISM. 
It can be seen that these results do span the entire range of observed [Sr/Ba] values, and 
even exceed the limits of the data. Obviously, in reality one would expect some mixing of 
shells to occur. As noted above, what is happening qualitatively is that the outermost shell 
is high in Sr and low in Ba, so produces a very large [Sr/Ba] value, whereas the innermost 
shells are low in Sr and high in Ba, so produce the lowest [Sr/Ba] values in the present 
GCE model. Additional smearing of these results might be expected from other effects that 
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are not included in our simple calculations, but the turbulence alone does appear to provide 
an explanation of the wide scatter observed in the data. 

If the atmospheres of these EMP stars are enriched by the tr-process, then the EOS 
dependence of the collapse time has a direct effect on the Sr and Ba production. Since the 
plotted values of [Sr/Ba] are already the maximum values, it would appear that in this 
model a very stiff EOS cannot produce the large observed values of [Sr/Ba]; a softer EOS 
is necessary. It's also noted that the upper limit of observed [Sr/Ba] values provides an 
observationally imposed lower limit on the softness of the EOS. 



5. Conclusions 

Our detailed inspection of Sr and Ba abundances in Galactic metal-poor stars indicates 
a clear cutoff in the distribution of [Sr/Ba] at [Fe/H]~ —3.6 and an upper bound as a 



function of [Fe/H] ([Sr/Ba] < -1.5-[Fe/H 



This work lends support to the suggestion (IBoyd et al. 



) as well as a low er bound at [Sr/Ba] ~ —0.5. 



2012) that the cutoff at [Fe/H] =-3.6 



for [Sr/Ba] >0.0 data for EMP stars can be explained by the tr-process, that is, by the 
collapses to black holes over a fairly wide mass range of progenitor stars. Furthermore, the 
results of our GCE calculations are found to predict the upper bound in [Sr/Ba] in these 
data. 

A fascinating aspect of this work is the possibility that observed minima in [Ba/Fe], 
[Sr/Fe], and [Sr/Ba] may be directly related to the stiffness of the nuclear EOS. The lower 
limit of [Ba/Fe] as a function of metallicity predicted in our model can be reproduced in a 
GCE model assuming a contribution to r-process elements from black holes produced in 
fallback supernovae. As shown in Figure [T]i, the softer the EOS, the lower the calculated 
[Ba/Fe] as a function of [Fe/H]. Thus the observed lower limit in [Ba/Fe] as a function of 
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metallicity appears to constrain the lower limit of the stiffness of the EOS. The tendency 



towards a soft 



( iFamiano et al. 



:qs, 



2006 



lowever, is consistent with 



Steiner. Lattimer. fc Brown 



jrior experiments and observations 



2010l ). While astronomical observations 



of neutron star masses are able to predict a lower limit of the stiffness of the EOS, this 
model suggests a method of determining the upper limit. 

Simulations also suggest that turbulence may explain the wide dispersion seen in these 
data sets, as the values of the predicted distributions do span more than the entire range of 
the data, and a more realistic model is necessary. 
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[Fe/H] [Fe/H] 



Fig. 1. — Results for (a) [Sr/Fe], (b) [Ba/Fe], and (c) [Sr/Ba] assuming a tr-process. The 
red dashed hne corresponds to a GCE model with no production in a tr-process, and the 
blue dashed line corresponds to a GCE model with a primary production from a tr-process 
for all stars with M>2OM0. Plot (d) shows GCE results for single-site tr-process production 
for turbulent ejection of specific shells assuming only those shells are ejected. 
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M/M 

sun 

Fig. 2. — Yields of Sr and Ba used for the GCE calculations in this model as a function 
of stellar mass and metallicity. The solid lines correspond to yields assumed for a primary 
r-process, while the dotted lines correspond to yields assuming a tr-process resulting from a 
collapse to a black hole assuming a LS220 EOS. 



